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Abstract 

Based on general symmetry principles we study an effective Lagrangian 
for a neutral system of condensed spin-1 deuteron nuclei and electrons, at 
greater-than-atomic but less-than-nuclear densities. We expect such matter 
to be present in thin layers within certain low-mass brown dwarfs. It may 
also be produced in future shock-wave-compression experiments as an effec- 
tive fuel for laser-induced nuclear fusion. We find a background solution of 
the effective theory describing a net spin zero condensate of deuterons with 
their spins aligned and anti-aligned in a certain spontaneously emerged pre- 
ferred direction. The spectrum of low energy collective excitations contains 
two spin-waves with linear dispersions - like in antiferromagnets - as well as 
gapped longitudinal and transverse modes related to the Meissner effect - like 
in superconductors. We show that counting of the Nambu-Goldstone modes of 
spontaneously broken internal and space-time symmetries obeys, in a nontriv- 
ial way, the rules of the Goldstone theorem for Lorentz non- invariant systems. 
We discuss thermodynamic properties of the condensate, and its potential 
manifestation in the low-mass brown dwarfs. 



1. Introduction and Summary 



The field theory method is an extremely powerful tool for studying properties of 
matter and radiation at various scales. The subject of the present work is deuterium 
matter at greater-than-atomic, but much less-than-nuclear densities. At such den- 
sities the neutral deuterium atoms are disintegrated into charged spin-1 deuteron 
nuclei and electrons, with an average inter-particle separation being smaller than 1 A 
but much greater than 1 fm. At temperatures of interest in this work (T <^ 10 5 K), 
the electrons form a degenerate Fermi gas. Deuterons, on the other hand, can be 
regarded as point-like spin-1 charged particles, which as argued below, will condense 
beneath certain temperature forming a quantum liquid. A field theory description 
of this condensed state of deuterons is the main topic of the present work. 

Recently, the condensation of charged spin-0 helium-4 nuclei at densities p ~ 
10 6 g/cm 3 was studied using the effective field theory method [1]. Such a dense 
system is expected to be present in helium-core white dwarf stars. It was shown [2] 
that these stars will cool faster because of the condensation of helium-4 nuclei (for 
a review see, [3]). 

Here we develop an effective field theory approach to describe the deuteron quan- 
tum liquid. This is interesting for a few reasons: (i) Densities at which deuterons 
could condense are much lower than those for helium-4 nuclei; as argued below, at 
p > 100 g/cm 3 , the deuteron condensate is expected to form upon gradual cooling 
of the system, (ii) The high-density deuteron matter is a promising fuel for laser- 
induced nuclear fusion. Hence, there has been a significant experimental interest 
in dense deuterium (as well as hydrogen) liquids (see [4] and references therein). 
Even though the present-day laboratory densities of shock-wave-compressed deu- 
terium matter are two orders of magnitude below 100 g/cm 3 , the experiments do 
see a transition to a liquid metallic state in which electrons get liberated from the 
atoms and the metallic conductivity is beginning to set in [4] . Given the importance 
of the subject one may hope the conditions necessary for reaching densities above 
100 g/cm 3 could be achieved in not too far future, (iii) The deuteron matter with 
p ~ 10 3 g/cm 3 is expected to be present in localized regions within low temperature 
(T <^ 10 5 K) low-mass brown dwarfs. These are compact astrophysical objects that 
have not managed - because of their low mass and temperature - to ignite nuclear 
fusion (i.e., they are "failed stars"). As we will argue, these low-mass, cold brown 
dwarfs should contain the deuteron quantum liquid in a thin layer in their interiors, 
(iv) Theoretically, the deuteron quantum liquid seems to be a richer system than 
its helium-4 counterpart. Due to deuterons being spin-1, the spectrum of small 
fluctuations of the deuteron liquid exhibits properties similar to those of antiferro- 
magnetic superconductors, as it will be shown in the paper. Moreover, the deuteron 
condensate represents an amusing example for which a consistency of counting of 
the Nambu-Goldstone modes in non-Lorentz-invariant theories [5] can be tested, (v) 
Our formalism can automatically be applied to condensates of other charged spin- 
1 particles, or to neutral spin-1 particles with contact, and/or spin-spin, and/or 
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magnetic dipole interactions, by taking various limits. 

Before we turn to the field theory description in the next section, it is instructive 
to summarize the conditions for the condensation to take place. 

At sufficiently high temperatures (but still significantly below the temperature 
of efficient thermo- and pycno- nuclear fusion), say at T ^ fO 5 K, the system of 
deuterons would form a classical Bose gas, even though the electrons may already 
be in a quantum degenerate state. At lower temperatures, however, transitions of 
deuterons to different states is possible - they may crystallize, or alternatively, may 
form a quantum liquid. 

Classical crystallization occurs at temperatures for which the Coulomb interac- 
tion energy of a pair of nuclei is about 180 times greater than the classical thermal 
energy of two degrees of freedom [6]. This leads to crystallization temperature for 
deuterons (see, [2] for recent discussions): 

T crys t — 1000 • p 1 / 3 K, (1) 

where density p is in units of g/cm 3 . At high-enough densities, however, before 
crystallization temperature is reached from above, the deuterons may undergo tran- 
sition into the charged condensate. This transition would take place around certain 
temperature T c for which the thermal de Broglie wavelength of the nuclei becomes 
greater than an average inter-particle separation between them. If T c significantly 
exceeds the crystallization temperature, then, before the system gets a chance to 
crystallize, the de Broglie wavelengths of the ions would begin to overlap leading 
to quantum-mechanical indistinguishability of the nuclei. If so, the standard ar- 
guments (see, e.g., [7]) suggest that the bosons would prefer to occupy one and 
the same quantum state, i.e., would prefer to condense. An estimate [1, 2] for the 
condensation temperature is 

T c ~ 350 • p 2/3 K. (2) 

The condensation temperature grows faster with p than the crystallization tem- 
perature does 1 . In particular, comparing (1) and (2), we conclude that for densi- 
ties greater than 100 g/cm 3 , temperature T c significantly exceeds T cryst , and hence 
deuteron condensation is favored upon gradual cooling of the system 2 . 

In the next section, based on universal symmetry principles, we formulate an 
order-parameter effective field theory describing the deuteron condensate. We dis- 
cuss the symmetry breaking patterns by the condensate, and count the broken gen- 
erators. In section 3 we obtain the solution that describes a spin-0 state of condensed 

1 For a discussion of uncertainties in (1) and (2), which should not be affecting much our con- 
clusions, see [8, 3]. 

2 At yet lower temperatures, when the crystalline state also becomes available, there may exist 
spontaneous transition from the condensate to the crystalline state. However, this process is 
strongly suppressed [8]. 
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spin-1 deuterons, and discuss small fluctuations about that solution. Interestingly, 
the fluctuations contain a massive mode, which is well below the ultraviolet (UV) 
cutoff of the effective theory, as well as two massless modes with linear dispersions. 
The massive mode is a reflection of the Meissner effect, while the massless modes are 
similar to antiferromagnetic spin-waves. Hence, the condensate has properties simi- 
lar to a antiferromagnetic superconductor, but has one order parameter. Moreover, 
the massless modes are Nambu-Goldstone (NG) bosons of spontaneously broken ro- 
tational symmetry. We discuss how counting of the NG modes agrees with the rules 
of Goldstone's theorem applied to non-Lorentz-invariant systems. In section 4 we 
study thermodynamic properties of a deuteron quantum liquid and briefly discuss 
its applications to low-mass and low-temperature brown dwarf interiors. Three ap- 
pendices summarize various important technical details which, if placed in the main 
text, would overwhelm the presentation. 



2. Effective Field Theory and Symmetry Breaking 

The condensed state of deuterons discussed in the previous section will be de- 
scribed here by a complex vector order parameter Wj(r, t), % — 1, 2, 3, that can slowly 

—1/3 

vary at scales much greater than the average inter-deuteron separation, d ~ J 
(hereafter J denotes the deuteron average number density, Jo ~ 10~ 6 MeV 3 ). A 
non-zero background value of the order parameter -its expectation value in the state 
with a macroscopic occupation (Wj(r, t)) - is a signature of the condensate. Small 
perturbations about the background, 5Wi = W,- L — (Wj), describe long- wavelength 
collective excitations of the system. 

In this section we will set up an effective Lagrangian for the order parameter 
Wi(r,t). Since the deuteron mass, ~ 1.876 GeV, is much greater than any 
energy/momentum scale that we'll be dealing with, the non-relativistic approxima- 
tion in ttih will be adopted (see Appendix A for the relativistic Lagrangian). The 
short-distance cutoff of the effective Lagrangian is set by the scale of interparticle 
separation d ~ J ^ 3 . As we will see, the effective Lagrangian will capture mass- 
less, as well as massive collective excitations of the system. The latter have masses 
~ giJo/n^H) 1 ^ 2 , where g denotes the gauge coupling constant. Because of the hier- 

1/3 

archy J ' <C to#, the above mass scale is parametrically supressed in comparison 
with the UV cutoff of the theory, d~ x ~ , and hence, can be calculated reliably 
within the effective Lagrangian approach. 

The effective Lagrangian has to be invariant w.r.t. all the internal and space-time 
symmetries appropriate for the physical system at hand. In particular, it should 
possess the local U(1)em symmetry of electromagnetism, the global subgroup of 
which guarantees the electric charge conservation. It should also posses a global 
U(l) s symmetry w.r.t. Wi(r,t) — > Wi(r,t)e t(3 with an arbitrary constant /3, as the 
latter guarantees conservation of the number of deuterons. The above two global 
symmetries, once put together, imply that the electron number is also conserved. 
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As to space-time symmetries, the Lagrangian should be invariant w.r.t. spatial 
rotations and reflections represented by the 0(3) group; it should also be Galilean 
invariant. Finally, in the limit when all the relevant coupling constants are set to 
zero, the equations of motion should reduce to the ordinary Schrodinger equation 
for the order parameter Wi. 

The Lagrangian that satisfies all the above conditions and describes the bosonic 
sector of the theory reads as follows [we use the "mostly-minus" metric convention] : 

A = -7 V" - kw* l D W t - W\D Q WiY) - ^—iD'Wj] 2 - 
4 2 2m H 

The covariant derivative is defined as = — igA^, g is a charge of W. The 
effective mass of a deuteron in the interacting system is denoted by m H . We will 
assume that this mass - being the highest scale - does not get renormalized too much 
by interactions and set m# ~ 2 GeV. The parameters k, X s , \ q are dimensionless 
coupling constants. The last term in the Lagrangian parametrizes the effects of the 
magnetic moment of charged deuteron nuclei for which we included the tree-level 
as well as loop-generated terms (the latter is parametrized by k). Furthermore, the 
A s - term describes the spin-spin interaction, while the X q - term parametrizes the 
contact repulsion of the deuterons 3 . 

We have included in (3) only leading order terms. There are higher dimensional 
terms in the effective Lagrangian that are consistent with the symmetries imposed. A 
systematic way to include those terms is to turn to the canonically normalized fields 
(replace W — > Wy/2m H ) ) and then write down the higher dimensional operators 
suppressed by appropriate powers of the short distance (UV) cutoff, d^ 1 ~ J^ 3 . 

1/3 

Because of the hierarchy J <C m#, the higher dimensional terms will not play a 
significant role for the low energy excitations of the system and will be ignored. 

In addition to (3) the full theory also includes terms describing electrons of some 
effective mass m e , for which we'll be using the Dirac Lagrangian 

£ 2 = + igA^ - mM- (4) 

The total Lagrangian reads as follows: C NR — L\ + £ 2 - For the reasons discussed 
above, we will also ignore higher dimensional terms which couple electrons directly 
to the order parameter (an example being W*Wipi(j). Even though we retain the 
Dirac Lagrangian in the canonical form (4), the electrons themselves are heavy and 
decoupled from the low-energy dynamics. What is not decoupled are the gapless 
excitations of the electrons near the Fermi surface. To begin with, we treat the 



3 Note that in a "fundamental theory" the k, X s , A 9 -terms need not be included in the tree-level 
Lagrangian as they'll be generated via loop corrections. In the effective Lagrangian, since we're 
not going to be calculating the loops, these terms were introduced with arbitrary coefficients in 
(3). 
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fermionic fluctuations in the Thomas-Fermi (TF) approximation, although this will 
be refined later on to include a loop correction due to those near-the-Fermi-surface 
gapless modes. In the TF approximation we average the fermionic four-current as 
follows: = (ip^tp) = ^o^o- This sets up a nonzero electron number- density J 
in the ground state. Note that this state does not break the global U(l) symmetry 
responsible for the fermion number conservation. 
As we will show in the next section the solution 

(A.H^C <wg = (o,o,o,t;), (5) 



where v = \J J /2m,H, represents a consistent classical background of the theory. 
It describes a state of condensed nuclei with zero net spin, in which the individual 
deuteron spins are aligned or anti-aligned with the preferred direction set by a vector 
v = (0,0, v). Such a spin alignment, and the resulting spin-waves discussed in the 
next section, are similar to those in antiferromagnetic materials. On the other 
hand, the properties of the electromagnetic field in this substance will be shown to 
be similar to those in superconductors. 

Note that a nonzero expectation value of Aq in (5) shifts the value of Fermi 
energy, however, this correction is typically of the order ~ A^IO -12 MeV, and will 
be ignored. 

Before turning to the next section, let us discuss the pattern of symmetry break- 
ing by the background (5) in a few specific cases: 

Case "a": A, ^ 0, A s 7^ 0, k ^ 1 . The internal part of the symmetry group G 
of the theory is given by Gj nt = U(l), whereas the relevant-to-us part of the space- 
time symmetry consists of the group of reflections and spatial rotations G st = 0(3) s 
(the Galilean and space-time translation subgroups of G are left unbroken by the 
condensate, and are ignored in the further discussions). Upon vector condensation 
(5), the symmetry group is broken down to spatial 0(2) s , corresponding to rotations 
and reflections in the plane perpendicular to v = (0, 0, v). Therefore, U (1) xO(3) s — > 
0(2) s , and the number of spontaneously broken continuous symmetry generators 
ends up being equal to 3. 

Case "b": \ q 7^ 0, \ s 7^ 0, k = 1 . The internal symmetry group gets enhanced by 
an additional factor of 0(3) j corresponding to a simultaneous internal rotation of Wi 
and W* by orthogonal three-by-three matrices. The diagonal part of 0(3), x 0(3) s 
gets broken by the vacuum (5) to the diagonal 0(2) d, whereas the orthogonal com- 
bination 0(3) a remains unbroken. This leads to the following pattern of symmetry 
breaking U(l) x 0(3)j x 0(3) s — > 0(3) a x 0(2) d , and consequently to three broken 
generators. 

Case "c": X q 7^ 0, X s = 0, k = 1 . Internal symmetry is enhanced further to 
Gi n t = 0(6). The symmetry breaking pattern is then 0(6) xO(3) s — > 0(5) xO(3);jn, 
where 0(3);jn is a linear combination of a subgroup of 0(6) and 0(3) s . This results 
in five broken generators. 



5 



The Lagrangian (3) describes a real vector field with two dynamical compo- 
nents, and a complex massive vector field Wi with six components three of which 
are heavy degrees of freedom (mass oc m H ) which should be ignored; this leaves al- 
together five potentially light degrees of freedom. Irrespective of the Case "a" , "b" 
or "c" described above, the U(l) electromagnetic symmetry is spontaneously broken 
by the condensate, hence, the longitudinal NG boson will become massive to serve 
as a helicity-0 component of a massive photon. Therefore, two NG bosons should 
remain and satisfy the counting rules for the different symmetry breaking patterns 
given in Cases "a" , "b" or "c" . At a first sight this seems unlikely, however, as we 
will see in the next section, that the dynamically obtained NG modes will obey the 
rules of counting of the Goldstone theorem for Lortentz-non-invariant models [5]; 
these rules are quite different from those in relativistic theories. 

Last but not least, in this and the next section we're discussing the effective 
Lagrangian at zero temperature (finite temperature effects in charged scalar con- 
densate have been recently studied in [9]), even though what we actually have in 
mind is temperature much lower than the deuteron condensation temperature, but 
higher than the room temperature. A caution has to be exercised in this regard - at 
a certain temperature T e <C 10~ 6 K (see the next section) there will be additional 
phase transition in which the Cooper pairs of near-the-Fermi-surface electrons will 
form via the Kohn-Luttinger mechanism [10] (see, [3] for a brief discussion of this 
effect in the presence of the BE condensate). Then, below T e a new NG boson 
will emerge as a result of spontaneously broken fermion-number symmetry trig- 
gered by the electron-electron order parameter. In conventional superconductivity 
an analogous NG Boson combines with a photon to form a massive photon. In the 
present case, however, this role is already played by the longitudinal mode associ- 
ated with spontaneously broken U(l) s boson-number symmetry discussed above 4 . 
On the other hand, since T e is close to zero for any system discussed in this work 
we'd expect temperature to be way above T e . Therefore, in our discussions we 
will not invoke the Cooper pairing mechanism with the fermion-number breaking 
condensate, and will ignore the finite temperature effects. 



3. Background and its Small Fluctuations 

Below we derive the background solution and study fluctuations about it. It 
is convenient to use a slightly different form of the Lagrangian (3) in which we 
"integrate in" an auxiliary field Wq, and also perform the following rescalings of the 
fields, W — > W^v^m//", A — > A/g. Then, the part of the Lagrangian relevant for 

4 In reality one should consider linear combinations of these two NG bosons; one of the combi- 
nations will become massive and another one will remain massless. However, since gJo/mn 3> T£, 
the mixing effects arc negligible. 
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studying the bosonic background and its fluctuations reads as follows: 

A = —F^F^ - ^E*^ + im H W*% - im H W % ^ + m 2 H W*W° - 

-J.N 1 - ^{W^Wif - ^-(e ijk W*W k ) 2 + iiAVW'F,,. (6) 
2 2 J 

where S^- = DiWj — DjWi denotes the gauged field strength of W, and the fermions 
have been treated in the mean- field approximation, as mentioned before 5 . There 
are no time derivatives of the W field in (6), and as noted before, it is an auxiliary 
variable which can be integrated out to get back the Lagrangian (3). 

Equations of motion for the fields Wq, W*, A , A* are straightforward to obtain 
from (6). They read respectively as follows: 

imjDWj + m 2 H W = 0, 
D%i - tm H D Wi + im H V i0 - (X q + X^W^W^Wi + X s (W j Wj)W* + 

iKW j Fij = 0, 

\diF j0 - 2m H W* i W i = J , 
J 2 ^F^ + iW*%i - %W 3 Y? 3l + m H (W*W + W t W*) - iK&{W*Wi - W*Wj) = 0. 
It is not difficult to check that these equations have the following classical solution 
A, - d,a = ^5°, W, = (0, e ia v), (7) 



where v is an arbitrarily-directed real vector with the norm v = \J Jq/ 2mu , and a 
is an overall phase of the vector field W^. As discussed in the previous section, we 
choose our background solution to be: 

A?; 2 

<i4 " ) = Wfr « = (0,0,0,*;). (8) 

This solution describes a state of a large occupation number of bosons, the net 
electric charge of which is compensated by the electron charge density Jo- The 
quantities that acquired the vacuum expectation values are gauge-invariant, as can 
be seen from (7). The gauge-choice given by (8), on the other hand, corresponds to 
the unitary gauge. 

To study the spectrum of small perturbations, it is useful to work in terms of 
the gauge-invariant variables 

A,-d^a 3 , \W 3 \, W 0)1 , 2 e- Ms , (9) 



5 A derivation of (6) as a non-relativistic limit of a relativistic theory in given in Appendix B. 
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where a 3 denotes the phase of W 3 (the same results would have been obtained 
by choosing the unitary gauge). As mentioned before, we treat the fermionic 
fluctuations in the Thomas-Fermi approximation: This can be accounted for in 
the equations of motion by increasing the coefficient of the A perturbation by 
TOq/ ' g 2 = (3n 2 Jq) 2 / 3 / 7r 2 in the ultra-relativistic regime for the electron Fermi-gas, or 
by m\j g 2 = m e (37i 2 Jq) 1 ^ 3 / tt 2 in the non-relativistic regime. 

The equations of motion for the excitations on the condensate, as well as the 
procedure of obtaining the linear spectrum are briefly summarized in Appendix C. 
The calculations are straightforward but very tedious. We find that the dispersion 
relations for the relevant, low-energy modes are given by the following: 

Longitudinal polarization of the massive photon. After symmetry breaking, one 
combination of the fields becomes the longitudinal mode of the massive photon, with 
the dispersion relation given by 

w? ~ 2v 2 g 2 1 + 2 I - . 10 
\ 8g z m z H v z J 

An identification of the corresponding mode is done by invoking the following argu- 
ments: first, we see that this dispersion is independent of the spin-spin interaction 
constant A s , as well as the magnetic moment interaction constant (1 — k), suggesting 
that this is a longitudinal mode. Second, in the limit of zero electromagnetic cou- 
pling, g — > 0, the dispersion relation for small k is linear in k and is determined by 
the "contact repulsion" constant X q , as it should be for a longitudinal mode 6 . Third, 
if in addition we take the limit X q — > we obtain the dispersion relation for a small 
fluctuation above a Bose-Einstein (BE) condensate of free particles, u± ~ k 2 /2m#. 

Transverse polarizations of the photon. A complete set of dispersion relations for 
the two transverse modes of the photon are given in Appendix C by 002,3 in (C-VII). 
In the limit when all the coupling constants are switched off we obtain ordinary 
photon dispersions. In the gauge theory the transverse modes acquire mass gaps, 
which can easily be seen by setting k — > 0. Expanding w 2 ,3 in (C-VII) in powers of 
the momentum for which |k| <C gv, one obtains 

u 2 2 ^2g 2 v 2 + k 2 + { -^£-^ (11) 

~ 2g 2 v 2 + k 2 + (1 ~* )2(k ' v) V (12) 
Amjj v 2 

Hence, the transverse modes are massive with the value of the mass coinciding with 
that of the longitudinal mode (10). A seeming independence of the last terms in 
the dispersion relations (11), (12) from the gauge coupling constant is due to the 
specific nature of the expansion |k| <C gv - the dependence on 1/ gv shows up in 



6 In the denominator of the second term in (10) we ignored terms proportional to m^k 2 and 
X q v 2 m,Q. These terms are always sub-dominant to the term retained in the denominator. 



8 



higher order terms. The exact expressions (C-VII) do no exhibit any pathological 
behavior. 

In the opposite limit |k| ^> gv, one may perform the l/m H expansion, leading 
to the following dispersions 

c 2 2 ^V + k 2 + 2sV%^k 2 , (13) 



I - v)-(k-- 
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ool ~ 2g 2 v 2 + k 2 + 2gV ^- ; v v " . (14) 

4m^ 

As we see, propagation of the transverse modes in general is affected by the magnetic 
moment interaction. However, cu 3 depends on the preferred direction of the spin- 
antispin alignment set by the vector v, while u 2 does not. As a result, for momenta 
perpendicular to the preferred direction, frequency u 3 is independent of (1 — re). 
On the other hand, when the direction of propagation coincides with the preferred 
direction, the expressions for 002 and W3 coincide too; this is a manifestation of the 
remaining rotational symmetry in the plane perpendicular to this direction. 

NG bosons of broken rotational symmetries. The two massless NG modes con- 
stitute the rest of the spectrum, their dispersion relations being given by 00^5 in 
(C-VII). In the zero coupling limit these reduce to the dispersion relations of small 
fluctuations over the BE condensate of free bosons, 00^ ~ /c 2 /2m#. In the inter- 
acting theory we perform the expansions similar to the ones used for the transverse 
modes, and obtain 

u 2 - ^ + A kV ^"^V (15) 
in the |k| <^ gv limit, while for the |k| 3> gv case we obtain 



^iir^S- 292 " 2 ^'- (l7) 



2 k 4 kV 2 2 (1 -K) 2 (k- V) 2 



4m^ s 2m 2 H 4m 2 H v 2 

A few comments are in order here: (i) As in the case of the transverse modes, 
seeming independence of the leading terms in NG dispersion relations (15), (16) 
from the gauge coupling constant is an artifact of the expansion in |k|/gt> 1. The 
spectrum is stable for re e [0,2] and A s > 0. (ii) The A s -dependent quadratic terms 
in the NG dispersion relations are due to the spin-spin interaction term (A s -term) 
in the original Lagrangian (6). (iii) The last, (1 — re)-dependent terms in (17) and 
(18), are due to magnetic moment interactions of oppositely aligned spins of charged 
deuterons. Since these interactions are attractive, the minus sign in front of the last 



9 



terms in (17) and (18) appears. In the exact expressions (C-VII) these negative 
terms can never overcome the positive ones, and the spectrum is stable, (iv) The 
dependence on the gauge coupling constant seems to disappear from the dispersions 
of the NG modes for k — 1. However, in the system of deuterons and electrons 
the effective constant of the spin-spin interactions, A s , and the anomalous magnetic 
moment of deuteron, k, will be proportional to positive powers of the gauge coupling, 
(v) As in the case of the transverse modes, for momenta aligned with the preferred 
direction, the dispersion relations for the two NG bosons coincide, as a reflection of 
the remaining 0(2) symmetry group. 

All the described properties above suggest that the obtained gapless states are 
NG bosons of spontaneously broken rotational symmetry. Their properties are very 
similar to spin-waves in antiferromagnets; in the present case they emerge due to 
spins of the charged nuclei. 

It is time now to check a consistency of counting of the NG bosons with Gold- 
stone's theorem for Lorentz non-invariant systems. This counting is rather different 
from that in relativistically invariant theories, and goes as follows [5]: if one de- 
notes the number of broken generators of continuous global symmetry by q, then 
the number of NG bosons should satisfy the following inequality 

q < m + 2n 2 , (19) 

where rii is the numbers of NG modes with frequencies proportional to an odd power 
of momentum in the long wavelength limit, (e.g., u oc |k|; these are called "type 
I"), while n,2 is the numbers of NG modes with frequencies proportional to an even 
power of momentum in that limit (e.g., u oc |k| 2 ; called "type II"). 

One can easily check, that (19) is remarkably satisfied by the NG modes found 
for the cases "a" , "b" , and "c" discussed at the end of the previous section. In the 
Case "a" or "b" , there are three broken generators one of which is Higgsed, hence, 
q = 2. This is matched by two type I NG modes, which saturate the inequality (19). 
If we impose a further condition that X s = - as in the Case "c"- then we end up 
with 5 broken generators one of which is Higgsed, leading to q — 4. On the other 
hand, both NG bosons found above end up being of type II in this limit, hence the 
inequality (19) is saturated again. 

One more interesting example is that of Case "a" with an additional condition 
X s = 0. The symmetry breaking pattern in this case is identical to that of Case "a" . 
However, the dispersion relations for the NG modes, (15) and (16), become of type 
II, i.e, we get rii = and n 2 = 2. Clearly, in this case the inequality (19) is not 
saturated, but is satisfied. 

As discussed at the end of the previous section, there will be one more longi- 
tudinal NG mode that emerges at extremely low temperatures (very close to the 
absolute zero) as a result of spontaneous breaking of the femion-number symmetry 
due to Cooper pairing of the near-the- Fermi-surface electrons via the Kohn-Luttinger 
mechanism. 
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We end this section by considering more refined effects of the near-the-Fermi- 
surface electrons. So far we have taken these effects into account in the Thomas- 
Fermi approximation. This approximation, however, fails to capture one of the most 
important properties of the fermion system, namely that the gapless modes near the 
Fermi surface can give rise to long range interactions for charged particles via the 
loop effects. 

One of the manifestations of this is the well-known Friedel potential [11], and 
the Kohn-Luttinger effect itself. It is important to see how these considerations- 
traditionally done for a system with fermions only - are affected by the presence of 
the spin-1 charge condensate (for spin-0 see [1]). 

To study these effects we consider a potential between two static charges (e.g., 
two impurities, such as helium-3 or helium-4 nuclei). The zero-zero component of 
the frequency-independent momentum-space Green's function in our case reads as 
follows (for simplicity, we take \ q = 0): 

k 2 

G(w = 0, k, A = 0) = — — r, (20) 

v ' ' > k 4 + 4M 4 ' v ; 

where M = (2g 2 m 2 H v 2 ) 1//4 . In the coordinate space this would give rise to an ex- 
ponentially suppressed oscillating potential (e _Mr /r)cos(Mr). We now include the 
contribution of the electron via one loop vacuum polarization diagram. This leads to 
two effects [11]: (i) the pole of the Green's function gets shifted (this is what leads to 
the Debye screening, and is the effect captured in the Thomas-Fermi approximation 
by the parameter m 2 ,); (ii) the gapless excitations near the Fermi surface, as in the 
case of the Friedel potential, lead to branch-cuts in the complex |k| -plane. These 
branch-cuts give rise to the long range piece in the position-space potential, which 
in the case with relativistic fermions (taken for simplicity), ends up being similar to 
the one obtained in [1] 

v nn ( e ~ Mr \ _l 4a em 4sin(2^r)\ 

Vstat = «cm<5iQ 2 cos(Mr) + . (21) 

\ r 7r M°r J 

Here denote the charges of static sources (impurities), a em is the fine structure 
constant and hp- Fermi momentum of the system. 

As in the case of the charged spin-0 condensate, the obtained long-range potential 
is suppressed in comparison with the standard Friedel potential. The reason for this 
suppression is the subtraction due to an exchange of the longitudinal polarization 
of the photon, as discussed in detail in [1]. 

The potential (21) is sign-indefinite and one may wonder whether its attractive 
part could trigger some instability near the Fermi surface to form Cooper pairs. In 
fact, the analogous oscillatory behavior in the interaction potential between electrons 
is what gives rise to the Kohn-Luttinger effect. However, as we see the potential 
is too small - a "typical" value for it is 10~ 16 MeV ~ 10~ 6 K, and therefore any 
tiny finite temperature will wash-out the effects of this potential. This small value 
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of the potential is also the reason for the temperature of the Kohn-Luttinger type 
transition to be very low - it should be that T e <C 1CT 6 K. 



4. Some Properties of the Condensate 

Already at densities p ~ 10 3 g/cm 3 (which is well above atomic densities for 
deuterium), the temperature of condensation significantly exceeds that of crystal- 
lization (the critical density, at which T c ~ T cryst is about an order of magnitude 
lower). Hence, according to the arguments given in section 1, at these densities the 
charged deuteron condensate is expected to form at T < T c ~ 3.5 ■ 10 4 K. 

In order to study thermodynamics of this substance at T T c , we need an 
expression for the temperature dependence of its specific heat {per ion). The latter 
will be determined by the massless NG modes 7 . The expression for the specific heat 
(in the units of the Boltzmann constant) is readily obtained from 

1 f kc 3 UJ 2 & ) I T 

Cv = j2^T J d * k 1)*T*' (22) 

Here k c = 2it/d represents the natural momentum cutoff in the effective theory. For 
densities mentioned above, k c ~ 10 _1 MeV. At the momenta close to the cutoff, 
the dominant contribution in the NG dispersion relation is given by the first terms 
in (17)-(18), which are of the order u 4 , 5 ~ 10~ 6 MeV ~ 10 4 K. 

We note that for the momenta greater than k = \/Al • 1CT 5 MeV, the k 2 /2m H 
term still dominates in the NG dispersion relations for most of the integration range 
in (22). A simple numerical analysis shows, that the contribution from the momen- 
tum domain k G (0, k ) is negligible as long as T ^> \ s v 2 /m H ~ A s • 10~ 3 K, which 
is assumed to be the case here. 

Therefore, the thermodynamics of the NG gas is determined by the free particle 
dispersion in the temperature range at hand, A s ■ 10 -3 K < T < 10 4 K. Using these 
estimates, we derive the expression for the specific heat per ion 

At temperatures T \ s v 2 /m H ~ \ s ■ 10™ 3 K, on the other hand, c v is determined 
by the \ s — term in the NG dispersion relations. Specific heat per deuteron in the 
ultra-low temperature limit therefore, is given by 

16tt 2 m 9 / 2 T 3 



15A^ /2 J 5/2 



(24) 



7 The gapped modes - two transverse and one longitudinal- are effectively "decoupled" due 
to their mass V2gv ~ 5 • 10 4 K being greater than the condensation temperature, while fermion 
contribution to the specific heat is negligible. 
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This is true as long as temperature is above the point T e where the Cooper pairing 
due to the Kohn-Luttinger (KL) effect will take place. Below this point there will 
be an additional NG boson contributing to the expression for specific heat. The 
latter will still scale as T 3 . Above or below the temperature of the KL transition, 
however, we expect that the heat transfer in deuteron condensate will be dominated 
by temperature- waves ("second sound") in the two-component phonon-superfluid 
medium. 

The densities of interest in the present work are expected to be present in certain 
astrophysical systems, namely in sub-stellar objects with masses low enough not to 
be able to fuse deuterium. An example is a low-temperature (T <^ 10 5 K), low-mass 
(M < 13 ■ Mj U pit cr ) brown dwarf. Such an object is supported by the electron Fermi 
degeneracy pressure working against gravity. An average interior density in the core 
could be close to ~ 10 3 g/cm 3 and the radius - of the order of Jupiter's radius 8 . The 
abundances of various nuclei (protons, deuterons and helium) in the core of such a 
dwarf are expected to trace the cosmological abundances of the corresponding chem- 
ical elements. In particular, the fractions of protons, deuterons and helium nuclei in 
the core should be similar to the fractions of the corresponding chemical elements 
in Jupiter's core N protons ~ 30 • N heHum ~ 10 5 • N deuteron . Furthermore, we assume 
a non-convective interior, and the chemical stratification of the elements (separa- 
tion of chemical elements) due to rotation of the object. Under these conditions 
the very center of the core would dominantly consist of helium sub-core of radius 
~ -Rjupiter/2, surrounded by a thin deuteron layer of ~ 2 km width, which is topped 
by a thick layer of protons, R pro ton ~ -Rjupiter/2. For temperatures around 10 4 K 
and densities at hand, the mean free path of the deuteron nuclei could be estimated 
to be ~ 30 m, suggesting that the smearing of the thin deuteron layer into the upper 
proton and lower helium layers is negligible. Also, the thermo- and pycno-nuclear 
fusion rates are small. Hence, we got conditions appropriate for condensation of 
the deuteron nuclei. The state of the core of this brown dwarf can be described 
as follows: in the sub-core the helium nuclei are crystallized; this is followed by a 
thin layer of a deuteron quantum liquid, which is topped by a gas of protons that is 
about to become quantum at slightly lower temperatures. Such an arrangement in 
the core should impact rotational properties of the brown dwarf. In particular, the 
rotation of the proton layer should not be expected in general to be in a phase with 
the rotation of the solid helium core, and therefore, certain irregularities in rotation 
might emerge. 
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Appendix A. Relativistic Description of the Condensate 

For a relativistic description of the system of deuterons in the background of elec- 
trons, we use an order parameter (an early treatment is given in [13]; for modern 
discussion, see [14]). The action should possess a local U(1)em, as well as Poincare 
invariances. The dynamics of the system is therefore described by the following 
effective Lagrangian 

Here = D P V V — D U V^ is the field strength (covariant derivative being defined 
as Dfj, = dfj, — iAfj), and J M represents the averaged fermionic current. It should 
be noted, that in this notation is rescaled by a factor of g, as compared to the 
physical normalization of the photon field. 

There exist three additional dimension-four, gauge-invariant operators which 
could contribute to the effective Lagrangian. These are the quartic self-couplings 
and magnetic dipole interactions of the form 

a = -\\i{v;v»y + l -\ 2 (y;v u y + mv;v v F^. (a-ii) 

For the purposes of studying the mass spectrum on the condensate, these terms are 
not of great importance and we will ignore them within the present section. They 
are however fully taken into account in the results of section 3. 
The equations of motion, following from (A-I) 

j 2 d^F, u + iiy^G^ - V»G%) = J v , D»G, U + m\V v = 0, (A-III) 

possess a solution 

A ll -d ll a = (m H ,Q), ^ = (0,e"V), (A-IV) 

where v is an arbitrary real vector with the norm v = a/ Jo/ 2m# and a is the overall 
phase of the vector V^. The norm of v is fixed by the first Maxwell equation that 
enforces the neutrality of the system in the bulk of the condensate. Our choice of 
vacuum corresponds to a = 0: 

(Afj) = (m H , 0), (V,) = (0, 0, 0, v) = v,. (A-V) 

This breaks the local U(l) symmetry, giving mass to the photon. It also breaks 
the spatial 0(3) group, down to 0(2). We should therefore expect two Nambu- 
Goldstone modes in the spectrum, corresponding to the latter symmetry breaking. 
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Hamiltonian of the system (excluding the contribution of the Fermi energy of 
electrons), on the above background 



n = -C = m H Jo, (A- VI) 

reduces to the rest energy of the vector bosons. 

Let us consider small perturbations of the background (A-V) in the equations of 
motion (A-III). By setting = v ^ + and = m H 5° + a M , one obtains 



□ + 2v 2 ) 77"" - — <W + 2u*V 



+[m H v v 5% - 2m H S^}(a^ + ^) = 0, (A- VII) 



[iv a d a rj^ - iv v d» - 2mv v 5% + m H 5 u ^}a^ + [(□ - 2im H d )v^ ~ d u d^ 

+im H (5 v Q d 11 + d u 5%) + m 2 H 5XH = °, (A- VIII) 

where r]^ = diag(l, —1,-1,-1) and all the derivatives act to the right. In order 
to obtain the mass gaps of perturbations on the condensate, we take the infinite 
wavelength limit (<9j — » 0, where % = 1,2,3). It is straightforward to show that 
the spectrum contains two massless modes, three masive photon degrees of freedom 
(mass V2gv), and three modes of mass 2m H which correspond to relatvistic effect 
of creation of boson-antiboson pair in the condesnate. 



Appendix B. Non-Relativistic Limit 

Restoring the factors of c, but still keeping h = 1, gives the following expression for 
the action 

S = J dtd 3 xC = J dtd 3 x ^—F^F^ - ^G; U G^ + m 2 H c 2 V;V» - 9 ~J^ ,(B-I) 

where the conventional normalization of the photon field is assumed, so that = 
8^ — igA^j c. The physical quantities are given by the following expressions 

E = --d t A-VA , B = VxA, p = — , f = J\ (B-II) 

c c 

where E and B are the electric and magnetic fields respectively, p denotes the 
number density of electrons and j is the corresponding current. The non-relativistic 
(magnetic) limit is taken by sending c to infinity, in a way that leaves physical 
quantities finite. This implies the finiteness of Aq, A/c and J°/c, and we neglect 
all quantities, suppressed by an extra factor of 1/c. It should be noted, that the 
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magnetic field happens to be enhanced by a factor of c. We switch to non-relativistic 
description by the following rescaling of the phase 

y» = e -im H cx f > w ^^ (B-III) 

where describes the non-relativistic vector field. The Lagrangian therefore be- 
comes 

1 X 

C = --F^F^ - + im H cW*% - irnncW 1 ^ + m 2 H c 2 W*W° - 9 -J^, 

where E M „ = D^W U — D u W fM . As a next step, we take the c — > oo limit keeping cW 
and Wj finite. This is reduces to rescaling W by a factor of c, and after setting 
c = 1 again, gives 

C NR = -\f^ - ^E*.E^ + im H W*% - im H W^ + m^ ¥° - 

-gJ^. (B-IV) 

Using the same reasoning, we could also include the rest of dimension-four operators 
(A-II) in the non-relativistic Lagrangian 

(C + C') NR = ~\f^ v - ^E*-E^ + im H W*% - im H W 1 ^ + m 2 H W*W° - 

-gJ^ - ^{W^Wif + 1a 2 (VU*W,) 2 + ////.il' -U •'/•;,. (B-V) 

One can integrate out the nondynamical W field from the last expression to obtain 
(3) 9 , which was our starting point for constructing the effective theory of section 2. 
The latter form of the Lagrangian makes the connection with Schrodinger's equation 
explicit. 

Appendix C. Spectrum of Perturbations 

In this appendix we briefly summarize the derivation of the full linear spectrum 
of perturbations on the deuteron condensate. The equations of motion for the 
perturbations {cro, ao, <J a , cr 3 , a a , a 3 } of the gauge-invariant variables {Wo, A , 
W a: I W 3 | , A a , A 3 }, which follow from the Lagrangian (6), read as follows: 

-im H d a (T a - im H d 3 a 3 - m H va 3 + m 2 H (r = 0, 

(C-I) 

X x 

— (-A + ml)a + —d (d a a a + d 3 a 3 ) + Am H va 3 = 0, 

g 2 g 2 

(Oil) 

im H d a a + [(-A - 2im H do)5 a p + d a dp] dp + d a d 3 a 3 - i(l - n)vd 3 a a - invd a a 3 + 

+\ s v 2 (a a -a* a ) = 0, 
(C-III) 

9 Parameters Ai and A 2 arc related to the counterparts of (3) by the following redefinitions 
A s = A2, A g = Ai — A2. 
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i 1 

-m H d 3 (a - (Jq) + (-A + 8$ + 2\ q v 2 )a 3 + -d 3 d a (a a + a* a ) - 2m H va = 0, 

(C-IV) 

- K)vd 3 ((r a - a* a ) - \d a d a + (— □ + 2v 2 )a a + \d a dpap + \d a d 3 a 3 = 0, 
9 9 9 9 

(C-V) 

m H v(a + o-q) - iKvd a (a a - <r* a ) ^d 3 d a + — (□ + <9|)a 3 + —d 3 d a a a = 0. 

(C-VI) 

Here a = 1, 2 and o" 3 is real. To obtain the spectrum of perturbations, we find the 
eigenvalues of the matrix, corresponding to the linear system (C-I)-(C-VI). After 
a tedious, but straightforward calculation, we find the following expressions for the 
spectrum 

2 2 2 (. k 4 + 2kVA n 

u>i = 2v g 1 + 



8g 2 m 2 H v 2 + (k 2 + 2v 2 \ q ) m\ 



1 



^2 2 4 = ^( k4 + 8g 2 v 2 m 2 H + 2k 2 (2m 2 H + v 2 \ s ) ± 
((Am 2 H (2g 2 v 2 + k 2 ) - 2k 2 v 2 \ s - k 4 ) 2 + 32^ 2 m> 2 k 4 (l - k) 2 ) 1 ' 2 ), 
uj 2 3 5 = ^V(k 4 + 8g 2 v 2 m 2 H + 2k 2 {2m 2 H + v 2 X s ) ± 



8m 



((Am 2 H (2g 2 v 2 + k 2 ) - 2k 2 v 2 \ s - k 4 ) 2 + 32^ 2 m^k 2 (k • v) 2 (l - /t) 2 ) 1/2 ). (C-VII) 

Here, the dispersion relation for the longitudinal polarization of the massive photon 
is given by uf; moreover, uj\ 3 correspond to the transverse polarizations of the 
photon, and are frequencies for two NG modes. The expressions for u 2 5 depend 
on the direction of propagation. Setting k G [0, 2] and \ s > 0, guarantees the 
absence of instabilities, as can easily be verified. 

It should also be noted, that the propagation of one of the NG bosons in the 
direction perpendicular to v is independent of the coupling constant g (putting aside 
the fact that \ s itself would be determined by g) 

^ 2 =iS; +A *S' for klv - (C ' VIII) 

This curious fact can be understood by a closer inspection of the system (C-VI). 
Integrating out the auxiliary variable <7o, as well as the nondynamical heavy modes 
cj) a = (<7 a + cr*)/2, one arrives at a simple equation for the light modes 

4m 2 <9 2 \ 

-A + 2\ s v 2 7r a = (1 - K)v(d 3 a a - d a a 3 ), (C-IX) 



where ir a = (a a — a* a )/2i. It can be easily seen for instance, that for k 2 = k 3 = 
(or &2 = d 3 = in position space) 7T2 represents the NG field, with the dispersion 
relation precisely given by (C-VIII), while 7Ti mixes with other bosonic fields of the 
theory, leading to the abovementioned asymmetry between the NG bosons. 



17 



References 



[1] G. Gabadadze and R. A. Rosen, JCAP 0902, 016 (2009) [arXiv:0811.4423 [hep- 
th]]; JCAP 0810, 030 (2008) [arXiv:0806.3692 [astro-ph]]; Phys. Lett. B 658, 
266 (2008) [arXiv:0706.2304 [hep-th]]. 

[2] G. Gabadadze and D. Pirtskhalava, JCAP 0905, 017 (2009) [arXiv:0904.4267 
[hep-th]]. 

[3] G. Gabadadze and R. A. Rosen, arXiv:0912.5270 [hep-ph]. 

[4] D. G. Hicks et al, Phys. Rev. B 79, 014112 (2009); W. J. Nellis, S. T. Weir, 
and A. C. Mitchell , Phys. Rev. B 59, 3434 (1999); L. B. Da Silva et al, Phys. 
Rev. Lett. 78, 483 (1997) 

[5] H. B. Nielsen and S. Chadha, Nucl. Phys. B 105, 445 (1976). 

[6] D.Q. Lamb and H.M. Van Horn, Ap.J. 200:306 (1975); S. Ichimaru, Rev. Mod. 
Phys. 54:1017 (1982); G.S. Stringfellow, H.E. DeWitt, W.I. Slattery, Phys. Lett. 
A 41, 1105 (1990). 

[7] A.J. Leggett, "Quantum Liquids," Oxford University Press (2006). 

[8] G. Gabadadze, Int. J. Mod. Phys. A 25, 627 (2010) [arXiv:0906.2406 [hep-th]]. 

[9] A. D. Dolgov, A. Lepidi and G. Piccinelli, JCAP 0902, 027 (2009) 
[arXiv:081 1.4406 [hep-th]]; Phys. Rev. D 80, 125009 (2009) [arXiv:0905.4422 
[hep-ph]]. 

[10] W. Kohn and J. M. Luttinger, Phys. Rev. Lett. 15, 524 (1965). 

[11] A. L. Fetter, J. D. Walecka, "Quantum Theory of Many-Particle Systems", 
McGraw-Hill (1971). 

[12] S.L. Shapiro and S. A. Teukolsky, "Black Holes, White Dwarfs, and Neutron 
Stars", John Wiley and Sons, (1983). 

[13] G. Velo and D. Zwanziger, Phys. Rev. 188, 2218 (1969). 

[14] M. Porrati and R. Rahman, Nucl. Phys. B 801, 174 (2008) [arXiv:0801.2581 
[hep-th]]. 



18 



